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Abstract. Due to the equivalence of the statistical ensembles thermostatic properties 
of physical systems with short-range interactions can be calculated in different 
ensembles leading to the same physics. In particular, the ensemble equivalence holds 
for systems that undergo a continuous phase transition in the infinite volume limit 
so that the properties of the transition can also be investigated in the microcanonical 
approach. Considering as example the spherical model the ensemble equivalence is 
explicitly demonstrated by calculating the critical properties in the microcanonical 
ensemble and comparing them to the well-known canonical results. 
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1. Introduction 

The properties of physical systems exhibiting a continuous (or discontinuous) phase 
transition are usually investigated in the canonical approach. However, for systems with 
short-range interactions the various physical ensembles such as the canonical ensemble 
and the microcanonical one are equivalent in the infinite volume limit so that the 
properties of continuous phase transitions can be investigated in the microcanonical 
ensemble as well.J The microcanonical description is based on the entropy as the 
thermodynamic potential and the physical properties of the system are deduced from its 
geometry. Microcanonical response functions, for example, are related to the curvature 
of the entropy surface. 

The description of phase transitions in the microcanonical formalism has gained 
growing interest in recent years ( [21 El IH El E] and references therein) . Apart from works 
about discontinuous phase transitions in microcanonical systems El El El HOI EEH El U2] 
second order phase transitions have been studied recently jU H31 EH UH E3 HH EH EH1- 
Ways to extract critical exponents from microcanonical quantities calculated for finite 
systems have also been suggested and applied successfully for various model systems 

\ It should be noted that the ensembles are inequivalcnt at a first order transition point even if the 
interactions are short-ranged. This subtle difference is exploited for the study of phase separation in 
microcanonical systems (e. g. pQ). 



Spherical Model in the microcanonical formalism 



2 



f3J EH ESB EH] • All these works basically concentrate on signatures of phase transitions 
in finite systems although some works investigated the general scaling behaviour of the 
entropy of the infinite system near a continuous transition [T^ l2~2~] . 

The main purpose of this paper is to demonstrate the equivalence of the 
microcanonical and the canonical ensemble in the infinite system for a concrete example 
of a model system with short-range interactions that undergoes a second order phase 
transition and can be tackled analytically. § To this end the properties of the spherical 
model at the phase transition point are investigated within the microcanonical ensemble. 
By doing so it is also shown exemplarily that the investigation of continuous phase 
transitions is not restricted to the canonical ensemble and can be carried out auspiciously 
in the microcanonical ensemble, too. There are only a few model systems with short- 
range interactions and a phase transition in the thermodynamic limit which, within 
the framework of the canonical ensemble, have been solved exactly. The zero-field 
Ising model in two dimensions and the spherical model in all dimensions belong to 
the class of models where the free energy density of the infinite system is known 
analytically. For all the author knows a microcanonical investigation of a model 
system with short-range forces exhibiting a non-trivial continuous phase transition 
in the thermodynamic limit has not been reported in the literature yet. Here an 
investigation of the interrelation of the canonical and microcanonical critical exponents 
of the spherical model is carried out. As the natural variables are different for the 
microcanonical and the canonical ensemble the critical exponents are in general not 
identical. The equivalence of the ensembles, however, leads to a relation between 
them. Note that this relation was discussed in for statistical ensembles whose 
thermodynamic potentials are connected to each other by Legendre transforms. The 
entropy as the thermodynamic potential of the microcanonical ensemble, however, is 
related, for instance, to the free energy by a Legendre transform and a subsequent 
partial inversion. The results for the microcanonical spherical model presented in this 
paper corroborate general considerations that have been based on scaling relations for 
the microcanonical entropy function [T3| I22j. 

The rest of the paper is organised as follows. In section |21 a brief introduction to the 
microcanonical analysis of physical properties of ferromagnetic systems is given. This 
section also sets up the notation and language used later on. The specific entropy of 
the ferromagnetic spherical model is calculated in section |3] in the macroscopic limit 
using the method of steepest descent. The critical properties of the spherical model 
are then analysed microcanonically in section 0] with special focus laid on the values of 
the microcanonical critical exponents. Section El contains some comments on the mean 
spherical model. The finding are summarised and compared to the canonical results in 
the final section |U1 

§ See for instance |231 1241 1121 12~5] for recent investigations of the question of the equivalence of the 
microcanonical and the canonical ensemble for systems with long-range forces. 
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2. Thermostatics in the microcanonical formalism 

The basic quantity in the investigation of the statistical properties of a finite magnetic 
system with N particles is the density of states 

n N (E, M) = J &T N 5{E - H{a))8{M - M{a)). (1) 

The Hamiltonian 7i provides the energy for any microstate a = (a*i, . . . , o~n) from 
the phase space of all possible configurations of the iV-particle system. The 
magnetisation M of the configuration a is measured by the operator Ai. Usually, the 
density of states is Laplace-transformed to yield the canonical partition function which 
determines the free energy. Note that the free energy is the thermodynamic potential 
of the canonical ensemble. In this paper, however, the density of states or equivalently 
the microcanonical entropy is directly analysed to deduce the physical properties of the 
spherical model. 

The microcanonical entropy density of a finite magnetic system with N particles is 
obtained from the density of states by taking the logarithm 

s N {e, m) = ^\ntt N (eN, mN) (2) 

where the energy density is defined by e = E/N, analogously, the magnetisation density 
is given by m = M/N. Note that units with &b = 1 are used in this work. The 
thermodynamic properties of the system in the thermodynamic limit are calculated 
from the entropy 

s(e,m) = lim s^(e,m). (3) 

The thermostatics of an infinite system can be investigated by studying the free 
energy as the thermodynamic potential with the temperature being one of the natural 
variables or alternatively by considering the entropy where the energy shows up as a 
natural variable [27| • In the following it is briefly summarised how the physical properties 
of statistical systems can be deduced from the entropy function s(e,m). The inverse 
temperature /3(e,m) and the magnetic field h(e,m) are basically given by the first 
derivatives of the entropy function: 

d 

P(e,m) = —s(e,m) (4) 

and 

Q 

h(e, m) = -(/9(e, m))" 1 — s(e, m). (5) 

am 

Note that the inverse temperature f3 and the magnetic field h show up as conjugate 
variables to the natural variables e and m in the microcanonical ensemble. 

The zero-field macrostate (e,m sp (e)) is defined to be the state with zero magnetic 
field h(e, m sp (e)) = for given energy e. The associated magnetisation m sp (e) is called 
spontaneous magnetisation. The corresponding inverse temperature is defined to be 
Po(e) = P(e,m sp (e)). The response functions of the system are related to second order 
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derivatives of the entropy function El E] The specific heat for the zero-field 
macrostate, for instance, is explicitly given by 



(A)(e)) 



2 



Co(e) = ~^7?7^- (6) 

^Po(e) 

In the vicinity of a critical point, which shows up at a critical energy e c , the physical 
quantities display power-law behaviour. Denoting the deviation of the energy e from the 
critical value e c by general physical quantity a has the form a(e) ~ |e| * E 

for small e. The singularity of a is characterised by the critical exponent x e in the 
microcanonical formalism. The critical exponent x t which characterises the singularity 
of a in the canonical formalism (i.e. a(T) ~ |T — Tel - ** with T being the temperature 
and T c being the critical temperature) is related to the exponent x e by the equation 

= Z • (7) 

1 - OC t 

The critical exponent of the specific heat in the canonical formalism is denoted by a t 
and has to be in the interval ]0, 1[ for relation (J7J) to be valid [I]. If the specific heat 
in the canonical formalism has a jump singularity, a logarithmic singularity or a cusp 
singularity at the transition point the microcanonical and canonical critical exponents 
are identical. 

As this paper considers the spherical model, which has a cusp singularity in the 
specific heat, the case of a cusp singularity is briefly sketched. In the vicinity of the 
transition point the canonical specific heat with a cusp at the transition temperature T c 
has the general form 

c (c) ~ A + B±\t\~ at (8) 

with t := T — T c and a negative canonical exponent a t . Integrating expression (|SJ) gives 

e^At ^-\t\ 1+ \^\. (9) 

i + kl 

The dominating term in the limit t — > is thus the linear term and one has e ~ t near the 
transition point. Therefore, physical quantities have the same qualitative dependence 
when expressed as functions of the reduced temperature t or the reduced energy e. 
The critical exponents are consequently identical for the canonical and microcanonical 
description. 

The exponent 5 e describes the relation between the critical magnetic field and the 
magnetisation at the transition energy (or temperature). The microcanonical exponent 
is therefore always identical to the canonical exponent St- 



3. Specific entropy of the spherical model 

In this section the density of states of the spherical model |2E] is calculated for finite 
systems. The spherical model exhibits a continuous phase transition and can be solved 
analytically for an arbitrary magnetic field in any dimension d. Therefore, its critical 
properties have been studied intensely within the canonical ensemble in the past (see 
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e. g. (213 EDI EH])- From the density of states one gets the specific entropy of the infinite 
lattice by carrying out the thermodynamic limit 0. The density of states will be 
evaluated for a hyper-cubic system in d dimensions with a linear extension L so that 
the system contains N = L d spins. The spin variables (7j £ R, % — 1, . . . , N, can take on 
any real value, but they have to satisfy the constraint 

JV 

5> 2 = iV - ( 10 ) 

8=1 

N 



The phase space of the spherical model is therefore the sphere of radius yN in 
The integration measure in the definition of the density of states is just given by the 
Lebesgue measure dlV = d N a. The Hamiltonian of the spherical model is given by 

H(a) = -jJ2^ (11) 

(id) 

with a positive exchange constant J, the magnetisation is 

M(o-) = J2°f (12) 

i 

The angular brackets indicate a summation over all neighbouring lattice sites i 

and j. The Hamiltonian (|11J1 together with the subsidiary condition (jlUJl defines the 
ferromagnetic spherical model with nearest neighbour interactions only. 
The density of states of the spherical model is generally given by 

tt N (E,M) = J d N a5(E -H(a))5(M - M(a))5(N ( i3 ) 

In view of the three delta functions the integral remains unchanged if one inserts the 
factors exp (aE — aH(a)), exp(bM — bM.(a)) and exp(ciV — c J2i °f ) with real a, b and 
c. Using the Fourier representation 

/dk 
— exp(ikx) (14) 

of the delta function one can rewrite the expression of the density of states as 

V N (E,M) =J^J^I^J d N ae W (A m:r (a)), (15) 

where the argument A of the exponential is given by 

Ap, g A<r) = (a + ip)E + (b + iq)M + (c + ir)N + Q Piq , r {a) ( i6 ) 

with the quadratic form 

Qp,q,r( a ) = ( a + W) 5^ 0i °i ~ ( C + ir ) ^ ~ ( b + Gi 

(id) « » 

= - a T Wa + v T a. (17) 

The last equality defines the matrix W, which describes the interaction of the spins, the 
vector v is the vector in which has N identical entries — (b + iq). The transposed 
vector of v e R N is denoted by v T . Introducing new spin variables \x G R w by 

^i = o--^W~ l v (18) 
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the form Q is given by 

Q P , q M = -fWn + ^v T W~ x v (19) 

in terms of the new variables \x (provided the inverse W~ x exists). The spin variables p, 
appear now quadratically so that one is left with a multiple Gaussian integral yielding 

fl N (E, M) =j^l^l ^exp($ p , 9 , r (E,M)) (20) 

with the argument 

I i 

$ Piq , r (E,M) = (c+ir)N+(a+ip)E+(b+iq)M+-v T W~ 1 v — lndet W. (21) 

Here the identity (det W)" 1 ^ 2 = exp(— \ lndet W) has been used. 

The interaction matrix W, that has been defined in expression (|17jh is a generalised 
cyclic matrix, which can be diagonalised by a Fourier transformation in any dimension 
d (e.g. [3U1 IHJ)- For the case of nearest-neighbour interactions only one gets the 
eigenvalues 

d 

A(y?i, . . . , ip d ) = (c + ir) - (a + ip)J ^ cos ipj (22) 

j'=i 

with (fj = (/ — 1)2tt/L, I = 1, . . . , L — 1 and j = 1, . . . , d. The logarithm of the 
determinant of W is now given by 

In det W N ~ 1 N J In Uc + ir) - (a + ip) J^cos^- J (23) 

[0,27r] d ^ j=1 ' 

where the summation has been replaced by an integral in the macroscopic limit of 
asymptotically large N (i. e. large L) . As v is an eigenvector of W with the eigenvalue 
(c + ir) — (a + ip)dJ it is also an eigenvector of the inverse W~ l with the eigenvalue 
((c + ir) — (a + ip)dJ)~ 1 . Therefore, one has 

-v^W^v = — — —. -r-rj,(b + iqf. (24) 

4 A((c + tr) - (a + ip)Jdy ' v ' 

Introducing the new integration variables z = a + ip, w = b + iq and u = c + ir the 
density of states is now given by 

a+ioo b+ioo e+ioo 

7v>i f dz f dw f du 

a—ico b—ioo 



n N (Ne,Nm)^ — — exp(iV0 w (e, m)) (25) 



with the argument 

<j>z,w,u(e, m) = ze + wm + u — 



1 1 (l^ ln (^-^E cos ^) 



w 2 



+ w^dTy (26) 
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This expression of the density of states for asymptotically large N can be used now 
to calculate the specific entropy of the infinite system by using the method of steepest 
descent |H2]- The argument of the logarithm in (j2Hj) must have a strictly positive real 
part for any possible value of the ipj. Therefore, the real part of u has to be positive. If 
Re z G ] — Rev,/ (dJ), Rev,/ (dJ)[ the function is analytic in z, w and u in this domain. 
Consider the function </> first for real values of z, w and u. For w — > ±00 one has <p ^ 00 
and for z — > ±u /(dJ) one also has — > 00. The function has therefore a minimum 
for a real wo, ^0 and vq with 2 G ] ~ v-o/ (dJ),vo/ (dJ)[. Due to the analyticity of (j> the 
integration path can now be deformed to have a = z , b = Wq and c = uq so that one 
obtains 

s(e,m)= lim — \nQ, N (eN,mN) = 4> ZOtWOtUo (e,m) (27) 

JV— >oo iv 

for the specific entropy in the thermodynamic limit by the method of steepest descent. 
The values zq, wq and vq as functions of the energy e and the magnetisation m are 
determined by the saddle point equations 

1 f d d ip JY^j C0S( Pj dJ 2 



dz 6+ 2 J (2n) d v- zJ^coscpj + A{v - zdJf W °' 

86 1 
ow 2(u — zdJ) 



and 



1 f d d <p 



dv 2 J (2ir) d v — zJ Y^j cos (fj A(u — zdJ)' 

4. Critical properties of the spherical model 



w 2 = 0. (30) 



In this section the physical properties of the spherical model are deduced directly from 
the specific entropy (|2*7|) . The main focus is laid on the possible appearance of a 
continuous phase transition signalled, for example, by diverging response functions. In 
particular, the character of the singular physical quantities at the phase transition point 
which shows up in dimensions larger than two will be scrutinised. As the entropy of the 
spherical model is an even function in m the discussion will consider only non-negative 
magnetisations. 

4-1. Discussion of the saddle point equations 

The microcanonical inverse temperature of the spherical model is given by 

3( e m ) = JL s ( e m ) = ^t + ^t^o d(j)_dwo d(f)_duo = d<p_ 
de ' de dz de dw de dv Q de de 

= z (e,m) (31) 

where the saddle point equations (|2*Sj). (|2*U|) and (|50jl have been used for the last but one 
equality. Similarly, the microcanonical magnetic field is determined by 

d 

h(e,m)f3(e,m) = —— — s(e,m) = —Wo(e,m). (32) 
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Defining the new function £(e, m) := u (e,m)/zo(e,m) the saddle point equation 
can be re-expressed as 

h(e, m) — 2m(C(e, m) — dJ). (33) 

Note that the variable £(e, m) has to be larger than the critical value ( c := dJ due to 
the restrictions for the argument of the logarithm in ()26|). The relation (|33J) can be used 
now to rewrite the other two saddle point equations (|2*K|) and (I37H) : 

B(em)\e + dJ m ^ = - 1 - ! 

m m ne + dJm\ 2 J (27r) d C(e,m) - J £ . cos ^ 

= : -P(C(e,m)) (34) 



and 



P(e, m) [I - m 2 ] = \ [ ^L- r ^ =: fl(C(e, m)). (35) 

2 J (27r) rf C(e,m) - J^.cos^- 

In the following the system for a fixed energy e is considered for the limit of 
vanishing (positive) magnetic field h — > 0+, i. e. for the zero-field macrostate (e, m sp (e)). 
The corresponding zero-field inverse temperature is defined to be (3o(e) = /3(e, m sp (e)), 
analogously Co(e) = m sp (e)). In the limit m —> m sp (e) (corresponding to h(e,m) — > 
0+ for fixed energy e) one gets the equation 

m sp (e)(( (e)-dJ)=0 (36) 

which has to be obeyed by the spontaneous magnetisation. This equation has the 
trivial solution m sp (e) = 0. A non-zero spontaneous magnetisation is only possible 
if the bracket in ()36|) vanishes, i. e. if Co = dJ. Therefore, an non-zero spontaneous 
magnetisation below some critical energy e c is only possible if Co(e) = dJ for e < e c . 

The possible emergence of a non-zero spontaneous magnetisation is further 
discussed in subsection 14.21 Before this consideration can be carried out it has to 
be investigated whether it is possible at all to have a critical energy e c larger than the 
ground state energy e g = —dJ so that Co( e ) = dJ holds below e c . To this end consider 
first the second relation (|33|) and assume that there exists a critical energy e c above 
which the spontaneous magnetisation m sp (e) vanishes. Then one must have Co( e ) > dJ 
for e > e c and the associated inverse temperature 0o(e) = ( 5(e,m sp (e) = 0) of the 
zero-field macrostates (e > e c ,0) is given by Po(e) = R(((e, 0)) (see relation (|33j) for 
the definition of R). In view of equation (|36|) a non-zero spontaneous magnetisation 
can appear for energies below the critical value only if Co(e) = dJ for all e < e c . At 
the critical energy e c one therefore has C( e c) = C and the associated critical inverse 
temperature f3 c is thus given by [3 C = -R(Cc)- In the limit ipj — > 0, j = 1, . . . , d, the 
sum £ . cos <fj tends to d and thus the vanishing denominator in (j33j) might cause a 
diverging integral for the limit Co — > Cc- This would have the consequence that the 
critical inverse temperature f3 c is infinite so that the model has not phase transition. 
To investigate this situation more explicitly consider the contributions to the integral 
which arise from small (pj. In the regime (fj — * the denominator can be approximated 
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Figure 1. Schematic depiction of the saddle point equation (|35|) for the zero- 
field macrostate. Dimensions d < 2 are illustrated on the left hand side whereas 
the right hand side shows the situation for dimensions d > 2. 



by !(y?i + . . . + (pj). Introducing polar coordinates for the d- dimensional 92-space and 
excluding a small sphere of radius S from the integration over <pj, j = 1, . . . , d, one gets 
the factor 



to whom the dominant small (fj contributions to the integral R are proportional. The 
modulus of the vector (</?i, . . . , ipd) is here denoted by ip. For one and two dimensions 
the limit ()37|) diverges and one has j3 c = 00 and the spherical model does not undergo 
a phase transition in these dimensions. From equation (|34j) it is evident that the 
associated critical energy e c = — dJ is the ground state energy of the spherical model. 
For dimensions larger than two the limit (|37p exists and hence the (critical) inverse 
temperature (3 C is finite and the model will have a phase transition at some energy 
e c > — dJ (see subsequent subsections). 

4-2. Spontaneous magnetisation 

This subsection discusses how a microcanonical spontaneous magnetisation can emerge 
in the spherical model for dimensions lager than two. Consider first the case of d < 2 
for the limit m — > m sp (e). The function R for the zero-field macrostate is shown 
schematically in figure ^ Below the critical value ( c = dJ the integral R is not defined. 
Above dJ the spontaneous magnetisation that is associated with £0 has to vanish in order 
to satisfy equation (|3*^1. On approaching dJ from above -R(Co) diverges and therefore 
one always has Pq = R(Co) (compare equation (J35j) ) and no spontaneous magnetisation 
can emerge. The situation for d > 2 is also displayed in figure ^ In contrast to the case 
d < 2 the integral R approaches a finite value R(( c ) in the limit — > dJ. If the inverse 
temperature (3q is chosen to be above -R(Cc) = Pc (or equivalently e < e c ) the saddle 
point equation ()35|) can only be satisfied if a non-zero spontaneous magnetisation m sp (e) 




(37) 



Spherical Model in the microcanonical formalism 



10 



shows up. || This is possible if Co( e ) — dJ for e < e c so that the subsidiary condition 
holds. Saddle point equation (j55]l now reduces to 



/ v / A)(e) - Pc (tiQ , 
m sp (e) = y-^- r . (38) 

The critical energy on the other hand is given by e c = -P(Cc) / Pc (see relation (f54"|) f° r the 
definition of P) and below e c the saddle point equation (|34|) for the zero-field macrostate 
is just Po(e + dJm1 v ) = e c /3 c . This equation can now be used to eliminate (3q from 
and one has 



m ^ e) = y^Tdi- (39) 

The spontaneous magnetisation is therefore characterised by the critical exponent 
/3 e = 1/2 for all dimensions d > 2. Note, however, that relation ()39|) is valid for all 
energies below the critical one. At the ground state energy e g = —dJ the spontaneous 
magnetisation is one as expected. 

At this stage it should be remarked that the entropy s(e,m) does exist for the 
spontaneous magnetisation m sp (e) although the argument of the logarithm in ()26|) then 
vanishes. The investigation of physical properties corresponding to derivatives of the 
entropy function gives sensible results for the limit m — > m sp (e). Note that similarly 
the free energy of the spherical model in the canonical formalism is also defined for 
zero magnetic field (i. e. for the spontaneous magnetisation) for temperatures below the 
critical one (e. g. [5T| 155]). 

4-3. Specific heat for dimensions d > 2 

The two saddle point equations (J54"|) and (|53]l for the zero-field macrostate contain the 
spontaneous magnetisation which can be zero or non-zero depending on whether the 
energy is above or below the critical value e c . Equation (|33j) can be used to eliminate 
the magnetisation from relation ()34|) yielding 



1 f d d (f dJ - J £\ cos (fj 



e = - =2 I dJ (40) 

2 J (27r)<<b(e)-J£i«*Vi/3b(e) ' 1 1 

Below e c one has Co = dJ so that 

e = \irr^- dJ = \ T ^ e )- dJ ( 41 ) 

where To(e) = l//3o is the actual temperature associated with the zero- field macrostate 
(e, m sp (e)). Thus, the microcanonical specific heat is constant, namely c = 1/2, for all 
energies below e c and all dimensions d > 2. 

|| In the microcanonical ensemble the energy e is the natural variable and the associated inverse 
temperature /?o( e ) for the zero-field macrostate is determined by the integral equations l|34|l and l|35|l . 
As the specific heat in the infinite system has a well defined sign, however, the inverse temperature /3o 
can be chosen first and from the saddle point equations the associated energy e can be calculated in 
principle, finally yielding the desired function /?o(e). 
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For arbitrary energies the microcanonical specific heat (JBJ) of the spherical model is 
given by 

1 f d d ip dJ - J cos tpj 



co(e) = \ J 



(27r)%(e)- JT.j<**<Pi 



+ -B (e) f dV dJ - J E 3 - cos p 3 - dCo 

+ 2 Po1 >J (27r)-(Co(e)-JE i cos^)>dA 1 J 

for the zero- field macrostate in view of equation ()4()jl . For e < e c the expression (|42jl 
reduces to the result already discussed above as d£o/d/3o = 0. In the limit ( — > Q+ from 
above the two integrals that appear in (|42|) both approach finite values which can be seen 
by a similar analysis carried out for the investigation of the behaviour of the integral R 
in this limit. The behaviour of the specific heat for the regime e — > e c + corresponding 
to Co —> Cc+ is thus determined by the behaviour of the derivative d£o/d/3o- 

As a first step for the analysis of the limit e — > e c +, i. e. e — > 0+, define for zero- 
field macrostate the new variables r := (3 C — /?o and £ := Co — Cc- Then r = R(( c ) — -R(Co) 
for energies above the critical value and one has the asymptotic relations 

( £( d - 2 )/ 2 if 2 < d < 4, 
r = J R(C c )- J R(Co)~ <^ -£m£ if d = 4, (43) 

^ if d > 4 

for the limit £ — > 0+ (2H1 EDI E^- The deviation of the energy from the critical energy 
is given by 

( = «.^S!.SVf + W)) - P,«). (44, 

The two differences -P(Co) — -P(Cc) an d R(( c ) ~ -R(Co)) however, have the same asymptotic 
behaviour for £ — > 0+ as the two integrals P and R differ only in the enumerator which 
does not alter the asymptotic behaviour for £ — > 0+. Thus, one has e ~ r near the 
critical point. Note that at this stage it is already obvious that the critical exponents 
of the spherical model are the same for both the canonical and the microcanonical 
ensemble. 

Using these results and the relation 

^Co = d£_ = 

d/3 d/? " dr 1 ' 
one obtains the following asymptotic relations for the microcanonical specific heat near 
the transition point: 

{1 /2 - A d ^ d -^/ 2 ~ 1/2 - if 2 < d < 4, 

l/2-4i/|ln£| -1/2-^/llnel if d = 4, (46) 

1/2 - fl d if d > 4 

where and Bd are some constants. In the analysis of the regime e — > e c the first 
integral in (}4*2*)) has been replaced by 1/2. Note that the correction term originating 
from this integral in (|42|) does not contribute to the leading asymptotic behaviour 
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of the microcanonical specific heat. For dimensions 2 < d < 4 the form ()46|) of 
the microcanonical specific heat is characterised by the negative critical exponent 
a E = (d — 4)/(d — 2) corresponding to a (right-sided) cusp singularity. For dimensions 
d > 4 the microcanonical specific heat has a discontinuity at the transition point. 

4-4- Susceptibility 

The susceptibility is generally defined by 

X = — =(—)~ 1 . (47) 
X dh \dm) 1 ' 

Focusing on energies e > e c and using the general relation l|H2j) between the magnetic 
field and the magnetisation one gets the microcanonical susceptibility 

UC) = 2Ue)~dJ) = k (48) 
for the zero-field macrostate. Expressed in terms of the energy deviation e one gets thus 
the asymptotic relations 

{ e ^d=2 if 2 < d < 4, 

(e/llne])- 1 if d = 4, (49) 
e" 1 if d > 4 

for the microcanonical susceptibility in the regime e — > 0+. The microcanonical 
exponent of the susceptibility in 2 < d < 4 space dimensions is therefore 7 e = 2/(d — 2) 
and in dimensions d > 4 one has the microcanonical exponent j e = l. 

4-5. Critical field 

At the critical energy e c the saddle point equation (|35|) reduces to 

/3(ec, m) = R(( c + h(e c , m)/(2m)) (50) 

where relation (p?2*j) has been used to eliminate ( from equation (|35|). Carrying out the 
asymptotic analysis of i? for the limit m — > m sp (e c ) = one gets the values 

*-{f «:>r 4 ' 

for the critical field exponent 5 e .^f 

Note that the analysis of the exponent S e is similar to the canonical analysis of St (see e. g. jHJ) an d 
is therefore not displayed here. 
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5. Mean spherical model 

The calculations presented in the above sections showed that the thermodynamic 
properties of the spherical model are equivalent in the canonical and microcanonical 
ensemble as far as the spontaneous magnetisation, the susceptibility, the specific 
heat and the critical field are considered. This is apparent at first sight from the 
corresponding expressions for these quantities. In the microcanonical ensemble the 
energy and the magnetisation are fixed to some values. These constraints are represented 
by the first two delta functions in the expression (|13[) for the density of states. These 
two constraints can be relaxed so that they are only satisfied on the average. This 
corresponds to the canonical treatment where the Lagrange parameters which are then 
introduced to satisfy the constraints on the energy and magnetisation turn out to 
be related to the temperature and the magnetic field, respectively. In a similar way 
the spherical constraint (jlOj) can be relaxed so that it is only satisfied on the average 
(the model is then often called mean spherical model). Lewis and Wannier used this 
relaxation to calculate the properties of the canonical spherical model [S3] . However, it 
turned out that the properties of the mean spherical model are only partially equivalent 
to the properties of the spherical model with the rigid constraint (fTUj) as quantities which 
are related to the fluctuations in the spherical constraint are different for the spherical 
and mean spherical model [SSI EH1 ESI EZJ • These differences have their origin in the fact 
that the Lagrange parameter that controls the constraint in the mean spherical model 
does not have a corresponding variable in the spherical model. Note that this is different 
for the Lagrange parameters temperature and magnetic field in the canonical ensemble 
which have corresponding variables, namely the energy and the magnetisation. 

A similar treatment can be carried out in the microcanonical ensemble. The 
spherical constraint can be relaxed and controlled by a Lagrange parameter so that it 
is satisfied on the average. Only two delta functions, namely those which fix the energy 
and the magnetisation, are left in the expression for the density of states and therefore 
only relations (|28|) and (|29|) show up as saddle point equations in the thermodynamic 
limit. However, the requirement that the spherical constraint is satisfied on the average 
leads to a subsidiary condition which fixes the introduced Lagrange parameter. This 
additional subsidiary condition is identical to relation (jBT)|) and therefore the investigated 
properties of the microcanonical spherical model and the microcanonical mean spherical 
model are equivalent. Nevertheless, the relaxation of the spherical constraint now allows 
for fluctuations in this constraint and having the canonical results in mind it might be 
expected that differences occur when comparing quantities that are related to those 
fluctuations. This interesting question goes beyond the scope of the present work and 
is left to future studies. 
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6. Summary 

The various statistical ensembles are equivalent for systems with short-range interactions 
undergoing a continuous phase transition in the infinite volume limit. Therefore, 
physical quantities can be calculated in different ensembles leading to the same 
thermodynamic properties. The properties of a model system undergoing a continuous 
phase transition in the thermodynamic limit can hence be directly deduced from 
the microcanonical entropy. As an example the spherical model was investigated 
microcanonically in this work to demonstrate these general properties and to provide 
an analytic treatment of a model system within the microcanonical formalism. The 
calculation of the microcanonical entropy as presented in section |3] somehow resembles 
the determination of the free energy in the canonical treatment (see e.g. [30]). This 
already hints at the equivalence of the microcanonical and the canonical treatment 
of the spherical model. However, in the microcanonial calculation further saddle point 
equations emerge due to the restrictions on the energy and the magnetisation represented 
by the delta functions in (|13|). To establish ensemble equivalence for the spherical model 
explicitly these additional saddle point equations have to be analysed as done in section 
□ 

The natural variables of the canonical and microcanonical formalism are different 
so that one gets in general different critical exponents for physical quantities. However, 
the canonical and microcanonical critical exponents are related to each other to ensure 
ensemble equivalence. In the case of a cusp singularity in the specific heat at the 
transition point the microcanonical and the canonical critical exponents are identical 
whereas the exponents are different for systems with an algebraically diverging specific 
heat (compare relation (jJJ)). The spherical model is a system with such a cusp singularity 
in the specific heat. The microcanonical critical exponents are therefore expected to be 
identical to the canonical critical exponents. This was indeed verified explicitly in this 
work. 

The microcanonical treatment of the spherical model on a infinite hyper-cubic 
lattice carried out in this work also demonstrates how properties of physical systems near 
a continuous phase transition point can be deduced directly from the density of states 
(or equivalently the microcanonical entropy) without going to the canonical ensemble. 
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